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A new variational technique is developed to investigate the polaronic fea- 
tures of the Holstein Molecular Crystal Model. It is based on a linear super- 
position of Bloch states that describe large and small polaron wave functions. 
It is shown that this method provides a very good description of the regime 
characterized by intermediate values of the electron-phonon coupling constant 
(the so-called intermediate polaron) for any value of the adiabatic parameter 
UJo/t. The polaron ground state energy in one and two dimensions is calcu- 
lated and successfully compared with the best estimates available providing a 
clear physical interpretation of the intermediate polaron. The band structure, 
the spectral weight of the ground state and the lattice displacement associ- 
ated to the polaron are also calculated and discussed. The new method has 
the advantage to require a very little computational effort. 
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I. INTRODUCTION 



Recently a large amount of experimental results, ranging from infrared spectroscopy to 
transport properties involving the colossal magneto- resistance and high superconductiv- 
ity, has pointed out the presence of polaronic carriers in doped cuprates and in the man- 
ganese oxide perovskites0'i. In particular several experiments have shown that in doped 
perovskite manganites Lai-x^xMnO^ {A = Sr, Co) there is a quite considerable coupling 
between the charge and the lattice degrees of freedom!. The lattice distortions associated 
with the Mn ions play an important role in determining the electronic and magnetic prop- 
erties of these compounds which have become the focus of the scientific interest after the 
discovery of the colossal magneto-resistance phenomena. Recently it has been found that 
in Lao. 756*00. 25Afn03 the metallic phase is characterized by homogeneously distributed in- 
termediate polarons while, above the transition from ferromagnetic metal to paramagnetic 
insulator, small and intermediate polarons coexist!. Also the measurements of the CuO 
distances in Laui^Sr^ i^CuOi crystal have pointed out, below lOO-ft', two conformations of 
the CuOq octahedra assigned to two different types of polaronJi. 

This large amount of experimental data has renewed the interest in studying the models 
of the electron-phonon coupled system and in particular the Holstein molecular crystal model 
that, for his relative simplicity, is the most considered model for the interaction of a single 
tight-binding electron coupled to an optical local phonon model. 

For the Holstein hamiltonian, beside the weak-coupling perturbative theoryi an ana- 
lytical approach is known for the strong coupling limit in the nonadiabatic regime (small 
polaron)!!. It is based on the Lang-Firsov canonical transformation and on expansion in 
powers of 1/A where A = Ep/zt is the dimensionless coupling constant, Ep, z and t be- 
ing, respectively, the small polaron binding energy, the coordination lattice number and the 
bare effective hopping integral (A represents the ratio between the small polaron binding 
energy and the energy gain of an itinerant electron on a rigid lattice). It is well known 
that both these analytical techniques fail to describe the region, of greatest physical inter- 
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est, characterized by intermediate couplings and by electronic and phononic energy scales 
not well separated. This regime has been analyzed in several works based on Monte Carlo 
simulationJ^il, numerical exact diagonalization of small clusterJi^, dynamical mean field 
theory0, density matrix renormalization group0 and variational approachesElQ. The gen- 
eral conclusion is that the ground state energy and the effective mass in the Holstein model 
are continuous functions of the electron-phonon coupling and that there is not phase transi- 
tion in this one-body systemlll. In particular when the interaction strength is greater than 
a critical value the ground state properties change significantly but without breaking the 
translational symmetry. 

Recently, results for the Holstein molecular crystal model have been presented by using 
the Global-Local variational methodEl. The comparison of the data obtained in one space 
dimension with the known approaches has shown that these results seem to provide the 
best estimates of the polaron ground state energy. They are highly accurate over a wide 
range of the polaron parameter space, from the non-adiabatic to the adiabatic, from weak 
to strong coupling limit. Nevertheless a solution of the Global-Local variational method 
for any particular k value {k is the wave number of the polaron Bloch state) is obtained 
by minimizing with respect to a very large number of parameters, that depends on the 
number of lattice sites and that increases dramatically with increasing the number of space 
dimensions from one to three. 

The aim of this work is to study the Holstein polaron features including ground state 
energy, polaron energy band, shape of the lattice distortion induced by the electron-phonon 
interaction and spectral weight of the coherent polaron band within a new variational ap- 
proach. It is based on two translationally invariant Bloch wave functions that provide a 
very good description of the two asymptotic regimes, the weak and strong coupling regimes. 
In this paper these wave functions are called large and small polaron. In the large polaron 
wave function the phonon distribution function takes into account the average effect of the 
correlation among the emission of successive virtual phonons by the electron and the spatial 
extension of the polaron is large compared with the lattice parameter of the crystal. In the 
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small polaron wave function the lattice polarization is confined to a region of the order of 
the unit cell and the polaron radius becomes of the order of the lattice constant. In this 
case all momenta of the Brillouin zone contribute to the polaron wave function and in the 
phonon distribution function the effect of the electron recoil due to the emission of virtual 
phonons is negligible. 

A careful inspection of these two wave functions points out that, far away from the two 
asymptotic regimes, they are not orthogonal and that the off-diagonal matrix elements of 
the Holstein hamiltonian are not zero. It is then straightforward to determine variationally 
the polaron ground state energy by considering as trial state the linear superposition of the 
large and small polaron wave functions. 

The comparison of our results with the Monte Carloi data and the ground state energies 
of the variational global local methodS shows that the proposed method provides a very 
good description of the polaron ground state energy for any value of the parameters of 
the Holstein model and confirms the existence of three regimes^: the weak coupling regime, 
characterized by polaron masses lightly heavier than the free electron mass and by dimension 
of the lattice polarization large compared with the lattice parameter; the strong coupling 
regime where the well-known polaronic band collapse takes place and the intermediate regime 
that is characterized by the crossover between the small and large polaron solutions. This 
regime is, therefore, well described by a wave function that is a linear superposition of Bloch 
states that describe the small and large polaron. 

We stress that the new variational approach provides a clear description of the Holstein 
polaron features in any regime and involves, for any particular k value, a very small number 
of variational parameters, that does not depend on the number of lattice sites. 

II. A NEW VARIATIONAL WAVE FUNCTION 
The model. The Holstein molecular crystal model is described by the Hamiltoniani: 
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H = H,i + Hph + Hj = -t J2 4cj + ^0 E (^h'T + E 4ci [Mge'^^-^^a^ + h.c. 

<i,j> q i,q 



(1) 



The units are such that h = 1. Since we will restrict ourselves to the single electron 
case we will not consider electron spin indices. The symbol <> in the first term of the sum 
in Eq.(P means that the summation is to be carried out only when i and j are nearest 
neighbours to each other. 

In the Eq.(|l]) c| denotes the electron creator operator at site i, the position vector of 
which is indicated by Ri, at represents the creation operator for phonon with wave number 
q, t is the transfer integral between nearest neighbor sites, ujq is the frequency of the optical 
local phonon mode and Mg indicates the electron-phonon matrix element. In the Holstein 
model (short range electron-phonon interaction) Mq assumes the form: 

Mq = -^u;,. (2) 



Here is the number of lattice sites. 
The small polaron 

When the value of g is sufficiently large the lattice polarization cannot follow the elec- 
tronic oscillations and, therefore, depends only on the average charge distribution of the 
electron. The wave function of the system can be factorized into a product of normalized 
variational functions \(p > and \f > depending on the electron and phonon coordinates 
respectively^!: 

I^(^) >= \v>\f> (3) 

where 

k>=E4lo>d0(i?m) (4) 

Rm 

and 1/ > has to be determined variationally. 

In the Eq.(^ |0 >ei is the electron vacuum state and (f){Rm) are variational parameters 
that satisfy the relation: 



Y.\<P{Rm)? = I . (5) 

The expectation value of the Hamiltonian ([l|) on the state (j^) gives: 

Rm,<S> 9 

(6) 

with 

p,^=M,Y.e''^-'^^\m)\' . (7) 

i 

In the Eq.(|^) the symbol <> in the summation means that 6 runs only over the nearest 
neighbours. 

The variational problem with respect to |/ > leads to the following lowest energy phonon 
state: 

1 „ —aff+h.c. , , ^ 

|/>=e^^L"o^ J|0>p^, (8) 
where |0 >ph is the phonon vacuum state, and to the following total energy: 

E^ = -t <P*{Rn.)<P{Rm-6)-Y^-^ ■ (9) 

Rm,<S> 1 

Even if this self-trapped state provides a good approximation of the ground state energy 
of the small polaron it is evident that the true eigenstate of the electron-lattice coupled 
system has translational symmetry. We construct translationally invariant Bloch states by 
taking a superposition of the localized states (0) centered on different lattice sites in the 
same manner in which one constructs a Bloch wave function from a linear combination of 
atomic orbitals. Then the trial wave function that accounts for the translational symmetry 
is given by (see also ref.l7): 

l4''^>=47Ee^'-'^"IV^?(^n)> (10) 

where 



Rn 



10 > 



ph 



(11) 



R„ 



and 



(12) 



This wave function is a sum of coherent states in the phonon coordinates, one for any 
particular lattice site. Since in a coherent state the emission of phonons occurs through a 
number of independent processes it is evident that in this trial state there is not correlation 
among the emission of successive virtual phonons. This physical assumption, i.e. that 
on every lattice site virtual phonons are emitted independently, is well-founded when g is 
sufficiently large but it is questionable for intermediate and small values of the electron- 
phonon interaction where the electron recoil kinetic energy plays an essential role. Moreover 
the wave function (|l^) does not contain states with real phonons. This indicates that the 
calculation of the polaron energy provides correct results only if the effective polaron band 
width A and the phonon energy uq satisfy the condition A < ujq. As it is well known, both 
these approximations limit the validity of the wave function ([T0|) to the strong coupling 
limit. 

The expectation value of the Hamiltonian ([l|) on the state ([lOD gives: 



Rn 



^ (j)*:{Rm)4>k{Rm - Rn - 5) 
Rm,<S> 



(13) 



J2rf:iRm)MRm-Rr, 



E 



Rn Rn 

2, , ^-iq-Rn _ ^iQ- {Rm- Rn) _ j^-^*^-iq-R„ 



(14) 



Moreover the calculation of the normalization factor < 1"^^*^ > gives: 



^4'*k{R7n)4'k{Rm - Rn) ■ 

Rm 



Rn 



(15) 



The next step is the determination of the variational parameters 0g(i?„). We note that 
if one neglects the spatial broadening of the electronic wave function, i.e. 0^(-R„) = q, 
the Lang Firsov approximation is recovered, i.e. the exact solution of the Holstein molecular 
crystal model for oj^/t oo. When the value of the adiabatic parameter uo/t decreases it 
becomes necessary to go beyond this approximation. In this paper we assume: 

Here /5g and 7^ are two variational parameters, is determined in such a way the Eq.(H) 
is satisfied and 5 and C, indicate, respectively, the nearest and the next nearest neighbours. 

This choice of the parameters and 7^ that takes into account the broadening of the 
electronic wave function in every lattice site to the nearest neighbours and to the next nearest 
neighbours allows to obtain a variational estimate of the ground state energy at = that 
is lower than the result of the second order of the perturbation theoryi: 

where Ep indicates the small polaron binding energy to the first order of the perturbation 
theory: 

^. = -E^ (18) 
and z is the nearest neighbour number. In appendix the calculation of the polaron band 



within the ansatz (|1^) for the variational parameters 0^(-R„) is reported. 

We end this section noting that this method can be systematically improved by adding 
further terms in Eg. ([T6|) . This allows to obtain better and better estimates of the polaron 
energy in the strong coupling limit. 

The large polaron 

It is well known that when the value of g is small the picture is quite different. As the 
electron moves through the crystal it exerts weak forces upon the ions which respond and 
move. This resultant ionic polarization will, in turn, modify the motion of the electron. 
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Then the particle must drag this polarization with it during its motion through the solid. 
This affects its effective mass that is weakly larger than that of a Bloch electronlll. In this 
weak-coupling regime is useful to adopt a variational approach similar to that of Lee, Low 
and Pines in the continuum approximation^. A possible choice for the trial wave function 



is: 



(19) 



Rn 



where 



>=4|0>eZ 



|0 >ph+T.d}ik)e-"^-''"a^^O 



(20) 



and 



Here Eiy[q) is the free electron band energy: 

d 

Eb{q) = -2t ^cos(gia) 

i=l 



(21) 



(22) 



where a is the lattice parameter and dq-{k) is a variational function that has to be determined 
by minimizing the expectation value of the Hamiltonian (|I|) on the state (0). 

\^^''\k) > has the right translational symmetry, i.e. it is a Bloch state with wave number 
k. This wave function represents an electron dressed by the virtual phonon field that de- 
scribes the ionic polarization. We note that the term in the square brackets of the Eq.(pOD 
allows a considerable advantage over the independent phonon approximation of Lee, Low 
and Pines. In the Lee, Low and Pines ansatz an important physical ingredient is missing: 
it does not take into account the fact that the polaron energy can approach ujq. On the 
contrary the wave function (|19|) contains this physical informationisl. In particular when the 
polaron excitation energy becomes equal to the energy of a longitudinal optical phonon, the 
band dispersion flattens and becomes horizontal. For these values of k the band has the 
bare phonon-like behaviour with very small spectral weight. For the same values of k and 



in particular at the edges of the Brillouin zone the main part of the spectral weight follows 
the bare electron bandi. 

The two wave functions (0) and (|1^), describing respectively the small and large polaron, 
differentiate mainly for the expression of the phonon distribution function. It is evident 
that, in spite of the assumption of no correlation (sum of coherent states in the phonon 
coordinate), the large polaron wave function takes into account the average effect of the 
correlation introduced by the electron recoil (Eq.plD, effect absent in the small polaron 
phonon distribution function. 

The expectation value of the Hamiltonian ([l|) on the state (p!9|) gives: 



<s> 
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12 



(23) 



q 



(24) 



where 
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Then dg is fixed by the condition: 



dE 



(0 



dd* 



Here E^^^ is the polaron energy in the weak-coupling limit: 



E 



(0 



<4 1-^14 > 

<^f\^pf> 



This procedure provides: 



dn 



y^, - % + 2^e-^^l'^'^'(^-™^^-'^) Eti cos (h - q^ 



(25) 



(26) 



(27) 



(28) 



10 



where 



i=l 



( ki \ ki 

cos Qi] a — cos — a 

\2 ^ 2 



ki 



sm Oi \ a — sm — a 

I 2 ^ / 2 



(29) 



and Xi{k), Zi{k), y{k) are variational parameters. 
Intermediate coupling 

For any particular value of t there is a value of the electron-phonon coupling constant 
(qc) where the ground state energies of the two previously discussed solutions become equal. 
Nevertheless the two solutions exhibit very different polaron features. In particular when 
the coupling constant is smaller than Qc the stable solution (the one with lowest energy) 
corresponds to the large polaron while for g > Qc it corresponds to the small polaron. 
Crossing the mass of the polaronic quasi-particle increases in a discontinuous way0. A 
more careful inspection shows that in this range of g values the wave functions describing the 
two solutions of large and small polaron are not orthogonal and have non zero off diagonal 
matrix elements. This suggests that the lowest state of the system is made of a mixture 
of the large and small polaron solutionsil. Then the idea is to use a variational method 
to determine the ground state energy of the hamiltonian (|1|) by considering as trial state a 
linear superposition of the wave functions describing the two types of previously discussed 
polarons: 



Al + Bl + 2A^B^S^, 

k: k. ft, ft, ft, 



(30) 



where 



wf > 



< 4" 1 4" > 



(31) 



and is the overlap factor of the two wave functions > and {ip^^'^ >: 



11 



In the Eg . (|30|) and are two additional variational parameters which provide the 
relative weight of the large and small polaron solutions in the ground state of the system 
for any particular value of k. 

In this paper we perform the minimization procedure in two steps. First the energies of 
the large and the small polaron wave functions are minimized, then these wave functions 
are used in the minimization procedure discussed in the present section. This way to pro- 
ceed simplifies significantly the computational effort and makes all calculations described 
accessible on a personal computer. 

It should be noted that the trial wave function ( PO] ) contains correlation between the 
emission of successive virtual phonons in the field around the electron since the phonon 
wave function is a linear superposition of coherent states for any particular lattice site. 
Then the wave function (^) recovers, in the weak and strong coupling limit respectively, 
the large and small polaron wave function, introduces correlation between the emission of 
successive virtual phonons by the electron and contains the important physical information 
that the quasi-particle becomes unstable when the polaron excitation energy equals the 
energy of a longitudinal optical phonon. 

The procedure of minimization of the quantity ^t^^jht^ with respect to Aj^ and Bj: gives 
for the polaron energy: 



E^, = , '2 ^^^^^ ^ (33) 



and 



S, 



^ ^ Eke - E^Sk 
% Er, - E? 



(34) 



k 

Here E^^ = (Ef + Ef) /2 and E^^ = '^f\H\4if > +/i.c.j /2. Finally the overlap 

factor and the matrix element of the hamiltonian between the two solutions {ip^^^ > and 
(g\ 

lipj^ > are, respectively: 
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:(<*f'i<'>)"'(<*f'i*r>)"' 

1 + 



and 



/2 



<s> 



-is) 



-Rr, 



-iq-Rr 



/2 



(35) 



(36) 



(37) 



III. NUMERICAL RESULTS 

In order to test the validity of our variational approach we recall the perturbative results 
both in the weak and strong coupling limits. From the weak coupling perturbative theory 
we geti: 



where 



\Mq? 



(38) 



S(A;, ikn) = ^ 

g- ikn -ujQ- Eb{k + q) 

while the second order perturbation theory in the strong coupling limit giveJ^: 



(39) 



/ 1 \ E 
E^r"^ ~ E„ ( 1 + — ) - 2te-»' cos k^ - 2— ^6"^' cos 2k^ 



(40) 
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in one dimension and 

~ (^1 + ^) - 2te-<^" (cos K + cos ky) 

E 2 

— 2— ^e~^ (cos2A;2: + cos2/c„ + 2cosA;2:CosA;„) (41) 

in two dimensions. 

In Fig.l. and Fig.2. we report the polaron ground state energy obtained within our 
approach (-E^^q in the Eg . (|33|) ) in one and two dimensions together with large and small 
polaron estimates (Eg. (^71) and Eg . (^6|) ) , on which our solution is based, and with the pertur- 
bative results. As it is clear from the plots, our variational proposal recovers the asymptotic 
perturbative results and improves significantly both variational estimates in the intermedi- 
ate region, where neither the perturbative methods nor the asymptotic variational ansatz 
give a satisfactory description. Moreover, our data in the intermediate region are in very 
good agreement with the results of two of the best methods available in the literature (see 
Fig. 3.): the Global Local variational method^! and the Quantum Monte Carlo calculationi. 
The agreement of our results with approaches numerically much more sophisticated indi- 
cates that the true wave function is very close to a superposition of the wave functions that 
we have classified as large and small polaron solutions. The very accurate choice of the 
variational wave function has allowed a dramatic simplification of the numerical problem. 

Within our approach we have also studied the polaron band both in one and two dimen- 
sions for different values of the electron-phonon coupling constant (Fig.4. and Fig.5.). As for 
the ground state energy our variational ansatz is able to recover all the properties expected. 
In the weak coupling regime, increasing the value of the wavenumber of the polaron Bloch 
state, Ej: increases until the excitation energy Ej: — E^^^ eguals ojq. When k is greater than 
this critical momentum the polaron becomes unstable to optical phonon emission and the 
dispersion curve bends over and becomes horizontal (this does not happen for t/uj^ < .25 in 
one dimension and t/uj^ < .125 in two dimensions). In the opposite regime the well-known 
polaronic band collapse takes place. Finally for intermediate values of the electron-phonon 
coupling constant the polaron band structure deviates significantly from both the dispersion 
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curves. In particular the strong coupling variational result underestimates the bandwidth 
and overestimates significantly the mass enhancement. 

From our results and in agreement with Romero et al.0 we find that there is not qualita- 
tive difference between the polaron features in one and two dimensions. In both cases, also in 
the adiabatic regime, there is a range of intermediate values of the electron-phonon coupling 
constant where a crossover takes place between the weak coupling regime, characterized by 
effective masses lightly heavier than the free electron mass, and the strong coupling regime 
in which the well known polaronic band collapse takes place. 

Another property of interest in studying the polaronic properties is the ground state 
spectral weight: 

Zu = \<i^k\A\^>? (42) 

where |0 > is the electronic vacuum state containing no phonons. is the renormalization 
coefficient of the one-electron Green function and gives the fraction of the bare electron state 
in the polaronic trial wave function. In Fig.6. and Fig. 7. we report the numerical results 
of Zg, at cuo/^ = 1) as a function of the electron-phonon coupling constant at = and 
as a function of the polaron Bloch state wavenumber for different values of g. In the weak 
coupling regime Zj:^^ is of order of the unity indicating that the polaronic quasi-particle is 
well-defined. The main part of the spectral weight is located at energies that correspond 
approximatively to the bare electronic levels. Instead at the edges of the Brillouin zone 
Z^ approaches zero. For these values of the wavenumber of the polaron Bloch state the 
main part of the spectral weight follows the bare electron band. Increasing the electron- 
phonon interaction Z^^^ decreases and approaches zero in the strong coupling regime. Here 
the carrier acquires large effective mass, the mean number of phonons in the cloud around 
the electron is very large and the most of spectral weight is located at the excited states, 
indicating that the coherent motion is suppressed rapidly with increasing the temperature^. 

Finally we consider the lattice displacement associated to the polaron formation. An 
estimate of the average deviation of the diatomic molecule on the site n + m from the 
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equilibrium position, when one electron is on the site n, is given by the function: 
where 

Sj;{Rm) = ^^^^^^^"'^^ . (44) 

Here M denotes the ionic mass and r^(i?„, R„i) represents the correlation function between 
the electronic density on the site n and the ionic displacement on the site n + m: 

Tj;{Rn: Rm) =< i>k\4iCn (a„+^ + a„+^) jV'g > . (45) 

In Fig. 8. we report the numerical results of the dimensionless quantity Sj^{Rm), at 
cuo/t — 0.5 and in one dimension, for different values of the electron-phonon coupling constant 

— * — * 

at A; = 0. In the weak-coupling regime Sj:^Q{Rm) decreases very slowly with increasing the 
value of m. This is consistent with the assertion that in this regime the extension of the 
polaron is large compared with the lattice parameter of the crystal. In the strong coupling 
regime Sj:^Q{Rm) is different from zero only for R^ — 0, i.e. the lattice displacement is 
different from zero only on the cell where there is the electron, indicating that the quasi- 
particle are extremely localized. Furthermore the crossover from large to small polaron is 
very smooth. 



IV. CONCLUSIONS 

In this paper a new variational approach has been developed to investigate the polaron 
features of the Holstein molecular crystal model. It has been found that a simple linear 
superposition of Bloch states that describe the small and large polaron solutions provides 
an estimate of the ground state energy that is in very good agreement with the best results 
available. It has been possible to identify a range of intermediate values of the electron- 
phonon coupling constant where a crossover takes place between the weak and strong cou- 
pling regime. Here the small and large polaron wave functions are not orthogonal and both 
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contribute to the formation of the so called intermediate polaron. We stress that the new 
variational approach does not require any significant computational effort to be implemented 
and involves, for any particular k value, a very small number of variational parameters, that 
does not depend on the number of lattice sites. 



V. APPENDIX 



In the Holstein model {Mg = uj^g/yN) the standard trigonometric integrals in the 
Eg ■ (p!3D , Eq.(0) and Eq.(|l5D can be performed analytically. In one dimension the polaron 
energy in the strong coupling limit assumes the following form [a = ujq = 1): 



E^^ = 7t (46) 



where 



< 4 l4 >= 1 + 2 COS A;, (2ag/3g + 2/?^7g) e ' ^ " " ^^^) 

/ o \ -g'^{a\-^fii+2'^i-2-ilal\ -g^ ( ai+2/3i+2'yi-2'yil3i] 

+2 COS 2A;^ + 27^a^j e V " "J + 4 cos 3fc^/?^7^e ^ " 

+2cos4fc^7fe V " h) ^ (47) 



, / ^9 \ -g^(ai+2l3i+2'ri-2all3i-2l3i'y^] 

-2t cos K i^oi^lk + + 1) e ^ ' ' 

/ \ -g'^(ai+l3i+2'^i-2-tlal\ 

-2t COS 2A;, (4/3^7^ + 2/3^ag)e ' \ >^ ' 

-2t cos (2a,^7g + /3| + 7|) e ' V ^^^^) 

-At cos 4A;^7g/5^e~^' {""t+^^t+^l) _ 2t cos 5fc^7?e~^' (^^E+^^E+^i) (43) 



and 



< ^f\H,n + >= -g' {4 + 2/5J + 27j 



2«|/5a: - 2/5f7fc - 2/5f a,^ - 27|/?g 
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2 cos krC 



In two dimensions we have 



2cos2/cre V * ^ 



< 1-0^ >= 1 + (cos A;^ + COS 

+ (cos + COS 2ky) (2/3| + 47I) e 
+ COS /ca; COS ky [^OL^^i + 8/3|j e V fc ^ ^ fc 



■4+4/34+474_2a|/32-4/3 



+ (cos 2/i;a; COS ky + cos 2A;j, cos k^ 8/?-^7-^e 
+ cos 2A;3; cos 



-t (cos /ca; + COS ky) (lO/3| + 47I + 2 + 8a^7gj e V u^'k ^fc 



-t (cos 2A;^ + cos 2ky) (4q;^/3^ + 16/3^7^) e"^' ("^^^s+'^d 



-t (cos 3A;^ + cos Sky) (2/3| + 47|) e ("s'^^^^+^^s) 
— t cos kx cos fcj^ (^16/5^0;^ 



^-,^(al+2Pl+4,l-2al,l) 



-t (cos 2A;^ cos ky + cos 2^;^^ cos k^) (l2/3| + I27I + 8a^-f^) e 



5 2/i;^ cos /Cj^ + cos 2ky cos /c, 
—t cos 2A;a; cos 2A;j^l6/3^7^( 
— i (cos 



3A;^cos A;^; + cosS/c^^cos /cj;) 8/3g7^e V* * 
—t (cos cos 2A;j^ + cos 3ky cos 2/ca.) 47^e V * * 



and 



< + («4 + 4^1 + 47I 



+9' 
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(cos kx + cos ky) e 



(2/5i + 47i) + 27|) - 4/3| - 87, 



(cos 2kx + cos 2ky) e 



32 4+3/34+274 



cos fca; COS /CyC 



-^^(-1+2/34+474-2724) 



(l6/5|7| - 8/?|7g - 8/3g7|) (cos 2fc^ cos ky + cos 2A;j, cos k^) e 



cos 2kx cos 2/i;^e 



(.4+4/34+374) 



(52) 
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FIG. 1. The polaron ground state energy {E^_q in the Eq.(|33|)) in one dimension (thick soUd 
Une) is reported as a function of the electron-phonon couphng constant for different values of 
the adiabatic parameter ojQ/t. The data obtained within the approach discussed in this paper 
are compared with the results of strong £'('^^) (Eq.(^), dashed-dotted line) and weak coupling 
perturbation theory E^'^^ (Eq.p^), thin solid line) and strong (dotted line) and the weak (dashed 
line) coupling variational estimates E^^ (Eq.(^6|)) and E^V (Eq.(p7[)). The energies are given in 
units of ujQ. 
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FIG. 2. The polaron ground state energy {Ej^_q in the Eq.(^3|)) in two dimensions (thick sohd 
hne) is reported as a function of the electron-phonon couphng constant for different values of 
the adiabatic parameter uJo/t. The data obtained within the approach discussed in this paper 
are compared with the results of strong £'(*^) (Eq.(^), dotted-dashed line) and weak coupling 
perturbation theory E^'^^ (Eq.p^), thin solid line) and strong (dotted) and the weak (dashed line) 
coupling variational estimates E^^ (Eq.(^6|)) and E^^ (Eq.(p7D). The energies are given in units 
of Uq. 
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FIG. 3. The variational results obtained within the approach discussed in this paper (Eq.(p3|), 
solid line) are compared with the data of the Global Local variational methodEl (diamonds), kindly 
provided by A. H. Romero, in one dimensions (Fig. 3a.) and with the energies calculated with a 
Quantum Monte Carlo algorithmi (circles), kindly provided by P. E. Kornilovitch, in one and two 
dimensions (Fig.3a. and Fig. 3b.) at oJo/t = 1. The energies are given in units of ojq. 
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FIG. 4. The polaron band structure (solid line) in one dimension at LOo/t = 1 is reported for 
different values of the electron-phonon coupling constant and it is compared with the weak (dashed 
line) and strong (dotted line) coupling variational estimates, E^^ and . The energies and the 
momenta are given in units of loq and n/a respectively. 
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FIG. 5. The polaron band structure along the diagonal {kx = ky) of the lattice in two dimensions 
at a;o/t = 1 is reported for different values of the electron-phonon coupling constant and it is 
compared with the weak (dashed line) and strong (dotted line) coupling variational estimates, £'1'^ 
and Ei'\ The energies and the momenta are given in units of loq and vr/a respectively. 
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FIG. 6. The spectral weight of the polaronic ground state in one dimension as a function of the 
electron-phonon coupHng constant at k = (Fig.6a.) and as a function of the polaron Bloch state 
wavenumber (in units of tt/o) for different values of g (Fig.6b., Fig.6c. and Fig.6d.) at ujo/t = 1. 
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FIG. 7. The spectral weight of the polaronic ground state in two dimensions as a function of the 
electron-phonon couphng constant at k = (Fig. 7a.) and as a function of the polaron Bloch state 
wavenumber (in units of tt/o) for different values of g (Fig.7b., Fig.7c. and Fig.7d.) at ujo/t = 1. 
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FIG. 8. The dimensionless quantity Sj:{Rm), at coo/t = 0.5 and in one dimension, for different 
values of the electron-phonon couphng constant at = 0: g = 1 (circles), g = 2 (squares), g = 2.2 
(diamonds), g = 2.5 (triangles). 
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